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We investigate a non-destructive measurement technique to monitor Josephson-like oscillations 
between two spatially separated neutral atom Bose-Einstein condensates. One condensate is placed 
in an optical cavity, which is strongly driven by a coherent optical field. The cavity output field 
is monitored using a homodyne detection scheme. The cavity field is well detuned from an atomic 
resonance, and experiences a dispersive phase shift proportional to the number of atoms in the 
cavity. The detected current is modulated by the coherent tunneling oscillations of the condensate. 
Even when there is an equal number of atoms in each well initially, a phase is established by the 
measurement process and Josephson-like oscillations develop due to measurement back-action noise 
^\ ' alone. 
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I. INTRODUCTION 



The experimental observation of Bose-Einstcin condensation (BEC) in dilute systems of trapped neutral atoms 
|y-[|] has stimulated a large research program on Bose-Einstein condensation of dilute neutral atom gases in confining 
otcntials. One aspect of BECs that has attracted much theoretical work is the idea of phase. Several papers 
7|, [i0||ll| ] have discussed the role of measurements in establishing the phase, in the form of interference between 
two condensates or Josephson-like coherent tunneling between the condensates. The latter situation is discussed in 
this paper, where we investigate homodyne detection of the output of an optical cavity containing one condensate 
in a double well system. The measurement process induces tunneling oscillations though back-action noise and thus 
induces a phase difference between the two separated parts of the condensate. In a self-consistent manner, the 
tunneling imposes a phase modulation on the light field, which is detected in the homodyne current. 



II. CONDENSATE MODEL 

The model of the condensate used here, namely the BEC in a double well potential, has been presented in previous 
papers |12(| , and so we only present a brief overview of the system here. The potential is symmetric with barrier height 
and well separation chosen so that only two single-particle states are below the barrier separating the two wells. This 
enables a treatment of the many body problem with a two-mode approximation. The resulting model is sufficiently 
simple to enable an analytic solution to be found for the semiclassical equations, and to permit a tractable numerical 
comparison of the semiclassical description with the full quantum dynamics. 

Consider a dilute gas of atoms moving in the double well potential, V(r) with, 



V(r) = b 



d 

2b 



+ -mu 2 {y 2 + z 2 ), 



(1) 



where the inter-well coupling occurs along x, and io t is the trap frequency in the y — z plane. This potential has 
elliptic fixed points at ri = +<7ox, T2 = — 9o x i where q 2 , = d/2b, at which the linearised motion is harmonic with 
frequency uiq — (Ad/m) 1 ' 2 . We set io t = wo for simplicity. It is convenient to scale the length in units of the position 
uncertainty in a harmonic oscillator ground state, ro where r$ 
B = (hcuo/8)(qo/r ) 2 



yJh/2mtjjQ. The barrier height is then given by 
The many-body Hamiltonian describing an atomic BEC in a confining potential is |13| | 



H(t) = / d 6 r 



V *2 



E/o 



2m 



v^ f • vv> + vfti> + —$$<$$ 



(2) 



where m is the atomic mass, Uq — 47r?i a/m measures the strength of the two-body interaction, and a is the s-wave 
scattering length, i/)(r, t) and i/j ( r , t) are the Heisenberg picture field operators which annihilate and create atoms at 
position r, and normal ordering has been used. 

For a suitable choice of B, only two energy eigenstates lie beneath the barrier which enables a many body treatment 
in terms of only two single-particle states. For details we refer to 12|. We now define the state uq(t) as the normalised 
ground-state mode of the local potential V^ (r), around the bottom of each well, with energy Eq, and define the local 
mode solutions of the individual wells lii^r) = wo( r — r i,2)- These local modes are approximately orthogonal with 
a first order correction (e 1 ) to orthogonality given by the overlap between the modes of opposite wells. The energy 
eigenstates of the global double-well potential may then be approximated as the symmetric (+) and asymmetric (-) 
combinations 



«±0) ~ -^ M r ) ± «2(r)] 



(3) 



with corresponding eigenvalues E± = Eq ± 7?., and 



TZ 



d 3 rul{r)[V{r) - V {2 \v - ri)]ua(r). 



(4) 



The matrix element TZ, which is of order e , describes the coupling between the local modes. The tunneling frequency, 
f2, between the two minima is then given by the energy level splitting of these two lowest states: 



a 2 ^2- 
n = 2K/h = ujAe 2r l. (5) 

2r 

In the two-mode approximation we expand the field operators in terms of the local modes and introduce the Heisenberg 
picture annihilation and creation operators 

c 3 (t) = J ' d 3 ru*(r)^(r,t) (6) 

so that [cj, c k ] « 5jk- The validity of this expansion if ensured when the overlap is small: 

ft fi 

p- = — « 1. (7) 

The ratio of the separation of the minima of the global potential V(r) to the position uncertainty in the state Uo( r ) 
can be as small as 2qo/ro — 6 (as in the simulations presented in this paper), and this condition is still satisfied. The 
many-body Hamiltonian then reduces to the following two-mode approximation: 

H 2 (t) = E Q {c\c x + c\c 2 ) + —{c x c\ + c\c 2 ) 



Tin 



c ((clfc? + (4) 2 ci) , (8) 



where k — U /2hV e j:f, and V~A = J d 3 r\u (r)\ 4 is the effective mode volume of each well. 

The two-mode approximation is valid when many-body interactions produce only small modifications of the ground 
state properties of the individual potentials. This is true when 

^ N \Uo\ 

hu> = — — » -^. 9 

2mr 2 V eff 

Using V e ff fa 87r 3 ' 2 Tq for this case, we obtain the following condition on the number of atoms 

N«^. (10) 

\a\ 

The values used in our simulations, namely ro = 5 /im, a = 5 nm and N — 100 satisfy this criterion. Thus the 
two-mode approximation is valid for small number of atoms compared to current experiments with N = 10 3 — 10 6 . 
The first term in Eq. H may be removed by transforming to an interaction picture, resulting in the Hamiltonian, 

Hi = ^(ci4 + c\c 2 ) + h K ((c\) 2 cl + (4) 2 c|). (11) 

A full quantum analysis of the quantum dynamics resulting from the many body Hamiltonian Eq(0) is not tractable, 
however considerable insight can be gained within the two-mode approximation. In |12f| an angular momentum model 
was defined which is equivalent to the Hamiltonian Eq(pT|). Using the transformations 

J z = -{c\c2 + c\ Cl ) (12) 

Jx = -^{c\c 2 ~ c\ci) (13) 

j y = l -{c\ Cl - c\c 2 ) (14) 

and setting c\c\ -I- c 2 c 2 = N = N (as the total number is conserved), the Hamiltonian becomes 

h 2 = nnj z + 2h K J 2 . (15) 

Here we have neglected terms proportional to N and -/V 2 since they merely correspond to a shift in the energy scale. 
The Casimir invariant is 




(16) 

This is analogous to an angular momentum model with total angular momentum given by j = N/2. 

The angular momentum operators have a simple physical interpretation. The operator J z corresponds to the 
particle occupation number difference between the single-particle energy eigenstates. For example the maximal weight 
eigenstate \j, j) z corresponds to all the particles occupying the highest single particle energy eigenstate, i/^ (x). The 
operator J x gives the particle number difference between the localised states of each well. In fact the x-component of 
the position operator in the field representation is 

x = -Jj- J '-c- ( 17 ) 

Thus the maximal and minimal weight eigenstates of J x correspond to the localisation of all the particles in one well 
or the other. 

III. HOMODYNE DETECTION SCHEME 

Figure 111 illustrates the system under investigation here. One of the wells of the double well system is placed in one 
arm of an optical cavity. The cavity is strongly driven by a coherent field at the cavity frequency. We assume that on 
the time scale of tunneling oscillations, the cavity is heavily damped. The cavity field thus relaxes to the steady state 
on a much faster time scale than the BEC dynamics. This enables us to make an adiabatic elimination of the cavity 
dynamics. The cavity field is assumed to be far off resonance with respect to a dipole transition in the atomic species. 
The effect of the atoms is then entirely dispersive and shifts the phase of the cavity field by an amount proportional 
to the number of atoms in the cavity at any particular time. If the atomic number in the cavity oscillates, so will 
the phase shift. Thus any tunneling of the condensate will be manifest in a modulated phase shift of the optical field 
exiting the cavity. To detect this phase shift we consider a homodyne detection scheme. The light leaving the cavity 
is combined with the reference beam and allowed to fall on a photodetector, which records the photocurrent. If there 
is a difference in atom number of the two condensates, then coherent tunneling can occur and the homodyne current 
will be modulated at the tunneling frequency. 

Assuming that the incoming light is detuned from any atomic resonance, the interaction Hamiltonian density is 

Hi = ¥{x)[H cm - hng(x)a f a]y(x) (18) 

where a, a' are the cavity field operators, g(x) is the intensity mode function and /i = f^/4A, with Rabi frequency 
17^ and optical detuning A. ^/(cc) and ^(x) are the atomic many-body operators, and H cm describes the centre of 
mass motion. 

Introducing the condensate field operators c, <y and averaging over the optical mode function gives the interaction 
energy: 

Hi = —hxa* ac\ci 

N + t " 

= — h— xa a — hxa aJ x (19) 

where x is the interaction strength. If the optical mode has a beam waist w, then the interaction strength can be 
written 



V2n 



(20) 



v/2(r /«;) 2 + 1 

For N — 100 atoms, x > 10~ 3 should give detectable phase-shifts (0.1 radians), and should be experimentally feasible. 
For example, with r$ as above, w = 30/xto, A/2ir = 100MHz, saturation intensity I s — 17W/m 2 , optical frequency 
u}/2ir = 3.8 x 10 14 -ffz, atomic linewidth 7 a /27r = 10 7 Hz and incident power P — 6mW, in a cavity 10cm long, then 
X ~ 10 -2 . Larger values of x may then be achieved by reducing the detuning or the incident intensity. 

The cavity is assumed to be driven by a strong coherent field of strength e and strongly damped at the rate 7. 
Hence the master equation for the whole system is, with h = 1, 



Ptot = -i£l[Jz, Ptot] - i2k[J%, ptot] 

- i{& ^-)W a , Ptot] +ixi a ■aJx,ptot] 

- ie[a f + a, p tot ] + - (2op tot a t - a) ap to t ~ ptot^a) (21) 

where the initial cavity detuning S = iV%/2 is chosen to remove the A-dependent linear dispersion. The optical field 
may now be adiabatically eliminated from the master equation |14| : 

p = -iCl[J z , p] - i2n[J x , p] 
+ i X \a\ 2 {J x ,p}-D[J x ,[J x ,p}} (22) 

where the coherent amplitude is a — —2ie/"f and D = 8x 2 e 2 /7 3 - The double commutator represents a decoherence 
produced by photon number fluctuations in the optical fields. It is a quantum measurement back-action term consistent 
with the interpretation that the optical field makes a measurement on the condensate. In fact this last term destroys 
coherence in the eigenbasis of J x and thus should inhibit tunneling oscillations. This is indeed true for the ensemble 
of measured systems described by the master equation. However, as we show below, it is not true for a particular 
realisation of a single measurement run. The ensemble-averaged effect of the measurement can be seen in the operator 
moment equations (for k = 0): 

< j x > = -Cl< J v > (23) 

< j y > = xM 2 < J z > +Q < X > -D < J y > (24) 



< j z > = - X \a\ 2 <J y >-D<J z >. (25) 



The terms with coefficient x\ a \ 2 produce a precession around the x-axis, which tends to inhibit coherent tunneling. 
The effect of these terms can be negated by adding a linear ramp, or tilt, to the double well potential. The D terms 
cause a decay toward the origin, indicating decoherence, as expected. If the system is started in a number state with 
an equal number of atoms in each well, then no tunneling will occur at all - these moments remain identically zero. 

When the wells are tilted, so that the precession around the x-axis is suppressed, we can obtain equations for the 
second order moments: 

< JJ > = -Q < A > (26) 
<J2> = 0<A> +2D(< Jf > - < J 2 >) (27) 

< P z > = 2D(< J 2 y >-< Jl >) (28) 
< A > = 2Q(< Jl > - < J 2 >) - D < A > . (29) 

where A = J x J y + J y J x . 

Thus even when the system is started with an equal number of atoms in each well, the unconditional evolution of 
< J x > and < J 2 > exhibit oscillations initially. For long times, the amplitude of these oscillations decays due to D 

and the system approaches the fixed point < J x >=< J 2 >=< J 2 >. From this we see that the condensate has on 
average a definite initial phase (which is clearly seen in the second order moments, but not the first order moments), 
which is determined by the initial state. In our simulations the initial state was chosen to be a eigenstate of J x , which 
is not the only state which gives an equal number of atoms in each well. Presumably in a real experiment, the initial 
state would be 

\^{T)>=e- iT ^\j,Q> x (30) 

where r is a random variable uniformly distributed on the interval [0, 2ii\. This then implies a random initial phase. 
A common technique for dealing with master equations describing open systems in quantum optics is to numerically 
simulate stochastic realisations of quantum trajectories. This method has already been used by several authors 
investigating the effect of measurement on the relative phase of BECs [0^1 , but these differ from the approach used 
here in that we monitor the homodyne detection current. The resultant stochastic process is a diffusive evolution rather 
than the jump processes which occur in the direct detection of atoms or individual photons. The quantum trajectory 
method is a very appropriate one for the situation considered in this paper. We have one two-part condensate system 
continuously monitored by the homodyne detection scheme. If there is any phase difference between the two parts 



of the condensate it will be established in each run of the experiment. The quantum trajectory method enables us 
to simulate each run of an experiment. The master equation however corresponds to an average over many runs of 
the experiment and many homodyne current records. For this reason moments calculated directly from the master 
equation will show no evidence of quantum tunneling if there is no initial phase difference between the condensates. 
In contrast, as we show, a single run of the experiment can establish a self-consistent phase difference even if no phase 
difference in present initially. Such a 'measurement induced' phase difference is manifest in a measurement induced 
tunneling current. 

The conditional master equation (that is the evolution conditioned on the measurement result) for the optical field 



undergoing homodyne detection is [TL6 17 



%) =^[a]p c + V 7 d ^H[a}p c (31) 

dt J field dt 

where dW(t) is the infinitesimal Weiner increment. In this equation, p c is the density matrix that is conditioned on 
a particular realisation of the homodyne current up to time t. Wiseman's superoperators are defined as 

T>[a]p = apa* — —(a*ap + pa* a) (32) 

Ti[a]p = ap + pa* — tr(ap + pa')p. (33) 

The stochastic Shrodinger equation, which describes the conditional evolution of the system is 

d|# c (t) >= dt[-iH - lpa)a + I(t)a]\V c (t) > (34) 

where H — VtJ z + 2k J% + Hi and |^ c (i) > is the unnormalised ket describing the conditional state of the system. 
The measured current is I(t) = 7 < a + a^ > (t) + y/j£,{t), where the stochastic term £(£) has the correlations 

< ew > = (35) 

<e(tU(0> =5(t-t'). (36) 

Thus we can see how the system evolution is conditioned upon the measured current. 
Adiabatic elimination of the optical field, using eq(17) gives 

d|* c (t) > = dt[-iH- ^-Jl +I{t)J x ]\* c {t) > (37) 

32y 2 e 2 4ve dW 

lit) = 6 -Ar- < J x > c +^=-^-(t). (38) 

Hence the oscillations in the occupation number between the two wells can be determined from the measured current. 
It is helpful to consider the Shrodinger equation for the normalised ket, which does not explicitly mention the 
detection current: 

d|* c (t) > = [-iHdt ~ %^{J X - < J x > c fdt + ^={J X - < Jx >c)dW}\V c (t) > (39) 

T V7 5 

The terms in the equation due to the measurement depend on the quantity J x — < J x > c . This is minimal in 
semi-classical type trajectories for which < jj > c factorizes to < J x >;?. Thus it may be expected that for some range 
of values of x, the stochastic measurement terms would drive the system towards an oscillating trajectory for which 

^ ° x ^ c— ^ Jx -^ c' 

IV. SIMULATIONS 

The results of the simulations are shown in figures 2 to 5. Time is plotted along the x-axis in dimensionless units 
(normalised by the inverse of the tunneling frequency SI). The strengths of the atom-atom collisions and the atom-field 
interaction were controlled by varying the normalised variables R — k/H and \ — X e /\A^T^ ^ ne parameters stated 
previously give the range of measurement strengths used in the simulations (10 -4 < \ < 0.1) when the power of the 
optical field is varied from 0.06mW to 6mW. The mass of the particles is taken to be m = 1.5 x 10~ 25 kg. 



The measured current gives the conditional dynamics of the the system. However, in the conditional results shown, 
we plot < J x > c , which is proportional to the current without the noise (eq |38|). For clarity, the other moments, 
namely < J y > c and < J z > c , are not plotted in the figures. Except when R is very large, < J y > c follows < J x > c , 
but with a tt/2 phase difference, and < J z > c remains at or close to zero. 

The dynamics of the unmonitored system when started with all the condensate in one well has been discussed in 
previous work |12]. Basically when there are no atom-atom interactions (ie R — 0), < J x > oscillates from —N/2 to 
+N/2. When the interactions are present but only weak, tunneling still occurs, but the amplitude quickly collapses 
due to nonlinear dephasing. The collapse is followed some time later by small revivals. There is a critical strength of 
collisions (RN = 1) at which the tunneling is suppressed. Above this value of R, the condensate is trapped in the well 
in which it started, with only very small oscillations occuring in < J x >. 

We expect to see similar behaviour in the current of the monitored system (eq^). When R = 0, < J x > c oscillates 
as before, for weak atom-light coupling (ie %iV ~< 1). For stronger measurements, the resulting back action can be 
seen in the current. For long times, the amplitude of the tunneling oscillations starts to fluctuate and a slow phase 
change is evident. In the case when atom collisions are present, the effect of the measurements is to halt the collapse 
of the oscillations seen in the unmonitored system. The phase changes are also more pronounced. The effect of the 
critical value (RN = 1) is seen in the suppression of the oscillations in the current above this value. Figure 2 shows 
the evolution of < J x > c for values of R above and below the critical value. 

If the system is started with an equal number of atoms in each well, then we expect no coherent tunneling in 
the absence of any detection apparatus. However, the presence of the field effects a measurement on the condensate 
system. This should establish a phase, which can be detected by measuring the output current I(t). The simulations 
of eq(|3^) show an oscillation in the current and, for the optimum interaction strength, this can be established in a 
few tunneling periods, for a small number of atoms. The results for the case where there are no atom collision, ic 
R = 0, are shown in figure 3, for various measurement strengths x- The growth in oscillations occurs because, for 
large enough x, < Jx >\ is driven to match < Jj > c , which typically has large oscillations. 

If, as in figure 3(a) and (b), the interaction strength is too small (xN < 0.1), then generally the fluctuations are 
not large enough to drive full tunneling and the current suffers small, rather irregular oscillations. However even for 
large x, when the oscillations in the current are established, such as in figure 3(c), they are not guaranteed to stay 
large in amplitude. This is because < J z > c undergoes what appears to be a random walk which, because of the 
Casimir invariant, directly affects the amplitude of the oscillations in < jj > c and < J 2 > c . Consequently, since the 
measurement locks < J x > 2 C onto the orbit of < J x > c , this changes the amplitude of the oscillations in the current. 
Because of the random nature of the orbit of < J z > c , the tunneling oscillations in the current over a certain time 
frame in a "good" run may be large, but in another with the same parameters the oscillations may be small and 
irregular in amplitude. 

When the measurement is quite strong (xN ~ 10), as in figure 3(d), the tunneling oscillations appear to be quite 
irregular. However, a fourier transform of the current picks out the tunneling frequency Q very strongly, so the 
fluctuations are mainly in the amplitude of oscillations, not so much in the phase. When the measurement strength is 
very large (xN = 40), < J x > c still indicates a tunneling from one well to another, but the oscillations are no longer 
harmonic and appear quite random, both in frequency and amplitude. 

The equations for the unconditional dynamics (eqs p3[ffi| ) show a decay in the oscillations of < Jj > and < J? > 
for long times, which increases with the measurement strength x- Figure 4 shows the unconditional dynamics for two 
different measurement strengths. No such decay is seen in the individual trajectories, but rather there is a diffusion 
in the phase of the oscillation over long times which accounts for the decay in the mean evolution. This change in 
phase appears to be most rapid over the periods when the oscillations are small in amplitude and most likely to suffer 
random fluctuations. 

The presence of atom-atom interactions increases the phase diffusion, even for quite weak collisions (RN < 0.1). 
Figure 5 shows the evolution of < J x > c for various atom-atom interaction strengths R, above and below the critical 
value. The amplitude is also more irregular, and the fourier transform of the oscillations no longer shows a clear peak 
at the expected tunneling frequency, but a group of random spikes centered on the the tunneling frequency. In figure 
5(b), R — 0.005 which is close to the value of R calculated from the parameters given in previous sections. Above 
the critical value of RN = 1, < J x > c suffers small, very irregular oscillations around the origin (figure 5(c)). This is 
quite different to the behaviour of < J x > c when the condensate was initially placed entirely in one well (figure 2(b)), 
in which case the condensated was trapped in the well it started in and the critical value of R marked quite a sharp 
boundary (or bifurcation) between two different types of behaviour. In this case where the condensate is distributed 
equally between the wells, the condensate is trapped in neither well, but remains across both, and the change in 
behaviour as the critical point is crossed is more continuous. As R increases past the critical point, we see a decrease 
in the overall amplitude of the tunneling and in its regularity. However, the critical value of R is still quite meaningful 



in this case an indication of when the strength of the atom collisions significantly suppresses the tunneling. 

V. CONCLUDING REMARKS 

We have shown that this homodyne detection scheme, for an appropriate choice of measurement strength \, could 
well be suitable to detect the relative phase, in the form of Josephson-like tunneling, between two condensates in a 
double well potential. The dynamics of the measured current reflect the tunneling of the condensate as well as the 
self-trapping effect caused by atom collisions. It also demonstrates quite vividly how a measurement can establish a 
(relative) phase in a system which initially exhibits no phase information. 
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FIG. 1. Schematic representation of the homodyne detection scheme to monitor the tunneling between two spatially separated 
condensates. One part of the condensate is contained in an optical cavity. The light in the cavity is well detuned from the 
atomic resonance. The output light from the cavity is detected by balanced homodyne detection. 

FIG. 2. Evolution of < J x >c in the monitored system, for N — 100 atoms all initially in one well and y — 0.01. In (a), 
R = 0.005 and in (b), R = 0.02. 

FIG. 3. Evolution of < J x > c in the monitored system, for N — 100 atoms and R — 0. In (a), y = 0.0001, in (b), y = 0.001, 
(c), x = 0.01, and in (d), y = 0.1. 

FIG. 4. Unconditional evolution of second order moments, for N — 100 atoms and R — 0. In (a), y = 0.01, in (b), y = 0.1 

FIG. 5. Evolution of < J x > c in the monitored system, for N = 100 atoms, and y = 0.001. In (a), R — 0.0001, in (b), 
R = 0.005, and in (c) R = 0.02. 
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